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ABSTRACT 

An effective chiral Lagrangian in heavy- fermion formalism whose parameters are constrained 
by kaon-nucleon and kaon-nuclear interactions next to the leading order in chiral expansion 
is used to describe kaon condensation in dense "neutron star" matter. The critical density 
is found to be robust with respect to the parameters of the chiral Lagrangian and comes 
out to be p c ~ (3 — 4)p - Once kaon condensation sets in, the system is no longer composed 
of neutron matter but of nuclear matter. Possible consequences on stellar collapse with the 
formation of compact "nuclear stars" or light-mass black holes are pointed out. 
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Introduction 



Kaon condensation in dense nuclear and neutron matter predicted first by Kaplan 
and NelsonQ has recently been reformulated in chiral perturbation theory (ChPT) to the 
leading chiral order by Politzer and Wise || and by Brown et al Q with results indicating 
dramatic consequences in stellar collapse || |5| and neutron star cooling (3|, Hj. In all 
these publications, the condensation process is mainly driven by the attraction gained as 
the K N sigma term n associated with explicit breaking of chiral symmetry 

Zkn ~ -(fh + m s )(N\{uu + ss)\N) (1) 

is "rotated out" by nuclear density, where fh = (m u + md) /2 is the mean up- and down- 
quark mass, m s the s-quark mass and u, s are the relevant quark fields. The K N sigma term 
([]]) which depends on the s-quark content of the proton could in principle be determined 
from experiment just as the irN sigma term could be determined from irN scattering. 
Unfortunately neither is the strangeness content of the proton known nor are K N scattering 
data accurate enough to determine the quantity T,rn- If we assume as an extreme case 
that {ss)n = (N\ss\N) ~ 0, then we obtain T,kn ~ l-3m n at the tree order with the chiral 
Lagrangian we write down below. If on the other hand we take (ss)n/ {(uu + dd))^ ~ 0.2 
which some of current hadrons models predict, we get T,kn ~ 3.7m n . In the literature, 
this range of magnitude of the sigma term has led to what appears to be a fairly robust 
prediction for kaon condensation at a matter density three or four times nuclear matter 
density po (~ ^m^). We shall consider this to be the range of values to keep in mind in 
what follows, although generally favoring the value of 

T, K N ~ Q^nN ~ 2m„- (2) 

using the empirical value of the ttN sigma term. The choice of this value will be justified 
later. 

The problem with this prediction based on the scalar attraction due to (|l|) is that by 
itself such an attraction is not compatible with available data on kaon-nucleon and kaon- 
nuclear interactions at low energy. As pointed out by several authors ||, there seems to 
be no serious need for a remnant of the scalar attraction in S-wave K^N and if ^ -nuclear 
scattering amplitudes with the bulk of the data adequately accounted for by standard p- 
and w-meson-exchange mechanisms 0. 

The purpose of this paper is to show that kaon condensation does still occur at the 
same low density, (~ (3 — 4)/9o), with an effective chiral Lagrangian that is compatible with 
low-energy kaon-nuclear interactions. The critical density could be lowered even further if 
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one invokes the previously proposed in-medium scaling (which we shall call "BR scaling") 
in the chiral Lagrangians M. We shall do this using heavy-fermion chiral perturbation 



expansion [ 10 1 to next to the leading order. We shall consider a chiral Lagrangian consisting 
of the octet baryons B = B a T a and the octet pseudoscalar mesons it = n a T a in heavy- 
fermion formalism (HFF) |j~0[| , first in free space and then in medium. We will focus on S- 
wave kaon-baryon interactions, so the decuplet baryons which we expect to be unimportant 
for this partial wave will not be included in our consideration. Although vector mesons do 
not figure explicitly in our effective Lagrangian, their chiral symmetry and scaling properties 
in nuclear medium will be invoked as in ref.[|| to infer the scaling behavior of the relevant 
term in the Lagrangian. 

Chiral Perturbation Theory 

In arranging terms in a consistent chiral expansion for mesons and baryons, it proves 



to be convenient to use the Weinberg counting rule |11] according to which an amplitude 
involving En number of external nucleon lines and Ek number of external kaon lines can 
be characterized by Q v in the amplitude where Q is the characteristic small momentum 
scale involved in the process and 

v = 2 + 2L - ^E N + Ok + \*H ~ 2) (3) 

~ i 

where L is the number of loops, the sum over i goes over all vertices, di the number of 
derivatives that act on the ith vertex and ni the number of nucleon lines attached to the 
ith vertex. In the absence of external fields, chiral symmetry constrains^] 

Pi = di + -m - 2 > 0. (4) 

Applied to K N scattering, we see that the leading term in this counting is given by L = 
and Pi = which is satisfied by a vertex with di = 1 and rij = 2, so the amplitude has the 
index v = 1. At the next order, we can have L = and one Pi = 1 vertex with an index 
v = 2. We note that no loops contribute to this order. One- loop terms contribute at the 
next order, say, v = 3 together with tree graphs involving one vertex with di = 2. In this 
paper, we will limit ourselves up to v = 2 for kaon-nucleon scattering and hence no loops 
need be calculated. Applied to kaon-nuclear scattering and to kaon condensation, this pro- 
cedure effectively takes into account density-dependent one-loop contributions consistently 



3 In the presence of slowly varying electroweak fields, the constraint is Pi > — f. This is important in 
nuclear exchange currents Jt2| . 
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with chiral symmetry as discussed in Q. This chiral counting is embodied in the chiral 
Lagrangian first written down by Jenkins and Manohar which we shall use: 



C = Co + £, (5) 



where 



f 2 

C = z-TrdnUdtrt + rTrM (U + h.c. - 2) 

+ TvBiv DB + 2DTvBS^ l {A fl ,B} + 2FTrB[A fl , B], (6) 
d = aiYiB (£M£ + h.c.) B + a 2 TrBB (£M£ + h.c.) + a 3 TrBBTr {MU + h.c.) 

+ Ci TtBD 2 B + c 2 TtB(v ■ D) 2 B 

+ d 1 TrBA 2 B + d 2 TiB{v ■ A) 2 B + d 3 TrBBA 2 + d 4 TrBB(v ■ A) 2 
+ d 5 Tr5STrA 2 + d 6 Tr5STr(w • A) 2 + ^TrSA^Tr^^ 
+ d 8 TrB(v ■ A)Tr(v ■ A)B + d 9 TvBA^BA^ + dioTri^ • A)B(v ■ A) 
+ hTvB{v -D)(S- A)B + / 2 Tr£(,S ■ D)(v A)B + / 3 T rJ B[S Q , S^A^A 
+ ••• (7) 
where is the four- velocity of the heavy baryon (with v 2 = 1) and 

L>^B = ^B+fy^j, (8) 

^ = i(C^C f + ^^0, ^ = ^(^^-^^0 (9) 

with £ 2 = C7. Here 5^ is the spin operator 5 M = 575 [A 7 M ] constrained to w • S = 0, Oj, 
Cj, cZj and /j are parameters to be fixed later. Except for some minor differences, we use 
the notations of Jenkins and Manohar |l(J where the advantage of using this Lagrangian for 
chiral perturbation theory is clearly pointed out. This Lagrangian was recently shown to 
be suited for describing chiral properties of nuclear systems by Park et al[12| where one can 



also find detailed discussions on how to systematically compute higher chiral-order terms 
using this Lagrangian. 

We now illustrate the counting rule discussed above in terms of the Lagrangian (]5|). 
The leading order KN (say K + N ) scattering is described by Co- Specifically the v = 1 
S-wave KN amplitude is given by iTrB[V°,B] which can be written explicitly as 

C u=1 = ^ (3(i\T 7 iV)K dt K + {Nt 7 °N) ■ Kt dt K) (10) 

with N T = (p n), K T = (K+ K°) and K dt K = K dt K - K dt K . For KN scattering, 
due to G-parity, the isoscalar term changes sign. In terms of an effective Lagrangian that 
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contains vector mesons such as hidden gauge symmetry Lagrangian of ref.[[D|, the first 
term of (|To| ) can be identified as the to exchange and the second term as the p exchange 
between the kaon and the nucleon. Thus we can think of the leading-order contribution 
as vector-dominated. This is in agreement with the standard phenomenological meson- 
exchange picture |7|. Now the next chiral order amplitude comes from £ (|7j). For S-wave 
KN scattering, the terms with the coefficients c% and fi do not contribute, so the relevant 
part of the Lagrangian simplifies (for S-wave) to 

C v=2 = ^L(NN)KK + ^{NtN)-(KtK) 



+^(NN)d t Kd t K + ^(NtN) ■ {d t Kfd t K) ( 1 I ) 



where 



^kn = -(^a 1 +a2 + 2a 3 )m s , (12) 

C = (13) 

D = j(d 1 + d 2 + d 7 + d 8 ) + ^{d 3 + d 4 ) + (d 5 + d 6 ), (14) 

& = ^(d 1 + d 2 + d 7 + d 8 ). (15) 

As stated, the Lagrangian (1^,11) at tree order gives the leading amplitude with v = 1 and 
the next-to-leading amplitude with u = 2. Loops with (^) give v > 3 and loops involving (|7|) 
even higher. Thus follows the statement that to order v = 2 there are no loop contributions. 
This makes the calculation simpler than in rnr or KK interactions where the next-to-leading 
order corrections involve both higher derivative counter terms and one-loop terms. This 
counting result in HFF is manifestly simpler than the relativistic formulation of ChPT [14] 
where numerous counter terms - none of which can be determined from other processes - 
intervene. 

In the strategy of chiral perturbation theory, the constants that appear in (JlOj) and 



(11) are to be determined from experiments. As mentioned above, the KN sigma term 
T,kn could, in principle, be determined from low-energy KN scattering data or from the 
expression (|l]) once we know the strangeness content of the nucleon (iV|ss|JV). At the mo- 
ment the sigma term is not known. We will take the value quoted above for our discussion. 
As for the coefficients D and D' we anticipate an important contribution from the 1/ms 
correction due to baryon-antibaryon pair terms. These "l/m^" corrections can be readily 
evaluated in the formalism (see Park et al |i"2]|). The argument is sketched in the Appendix. 
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Here we simply give the results 
1 r. 



Dl « 

m 48 

6', « -i- r 

m 48 



D + 3F) 2 + 9(D - F) /m B rj -0.12/m B ~ -0.024fm, 
D + 3F) 2 - 3(D - F) 2 1 « -0.086/m B « -0.017 fin (16) 



where m B is the (centroid) baryon mass which we take ~ 1 GeV. We have used here the 
tree-level fits F = 0.44 and D = 0.81. There will also be contributions that arise from 
degrees of freedom whose mass scale is higher than the chiral expansion scale A x ~ 1 GeV. 
These can appear as a counter-term contribution and cannot be calculated in the scheme. 
They should be determined from experiments. In some cases such as in irir scattering, those 
constants can be saturated by resonances (for instance, the constants Lj in the Lagrangian 
of 0(Q 4 ) in the pion sector are dominated by the vector mesons p, a% etc.) but here while 
this is plausible, we have not succeeded in finding such a simple mechanism. 

On-Shell Constraints 

What can we say about the constants D and D' from experiments? 

For the purpose of determining these constants from experiments, it is simpler to 
look at K + N scattering. The K~N scattering is somewhat more delicate because of the 
resonance A(1405). It has been suggested g, 0, |] that the presently available data on K^- 
nucleon and .RT^ -nuclear scattering indicate that the bulk of the data can be understood 
reasonably well by the leading v = 1 term ( p!o| ) with the v = 2 terms effectively suppressed. 
Let us see what this means in the framework of chiral perturbation theory. At present 
the data on threshold K + N scattering, particularly the isoscalar amplitudes, are not good 
enough to allow a precise determination of the constants. We can nonetheless make a simple 
analysis which is still meaningful as described below. 



From the effective Lagrangian (10) and (|il|), we can immediately write down the 
expressions for scattering lengths 

°^ = a f2M T 1 -A-m K + V KN + C+(D + &)m 2 K ), (17) 

47T/ Z (1 + ttik /tub) v 1 

a ?=° = A f 2n / / S (XKN-3C+(D-3D')m 2 K ). (18) 

47r/ z (l + mx/mB) v ' 

where af=g \ are the S-wave scattering lengths for K + N scattering in isospin / = 0,1. 
Given experimental values for aj=oi and knowing T,kn, these equations could determine 
the constants D and D' . 
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Scattering lengths have recently been obtained from a rather complete analysis by 
Barnes and Swanson [15]. They find 



af^ « -0.31fm, 

afl « -0.09 fm. (19) 

Whereas afj is reliably determined, af= Q is highly uncertain: It can differ between and 
-0.2 fm |l|g. 



As long clS dj q IS small in magnitude compared with a|^, the nuclear scattering 

of kaons should be insensitive to afj^. This suggests a procedure which minimizes the 
dependence on the precise value of af= . The isospin averaged amplitude 

a = 3 -afJ + \an (20) 

will occur in kaonic interactions in nuclear matter. In the case of neutron star matter, 
another combination occurs, but we will show that this latter case is insensitive to how we 
handle the on-shell kaon scattering. From (19), we find the empirical value 

flemp ~ —0.255 fm, (21) 

which is rather insensitive to the precise value of aj=o within the given limits. 

Miiller-Groeling et al Q have analyzed KN scattering in the framework of boson 
exchange. They employ vector mesons with roughly the same SU(3) coupling coefficients 



as are implicitly incorporated in eq.(10). Whereas the scalar attraction between kaon and 
nucleon is built into the chiral Lagrangian ( |Tl|) through the T,kn term, Miiller-Groeling et 
al H include an exchange of an explicit scalar particle with 

^tnn^kk^tt m 0.9 (22) 



and m a ~ 600 MeV. Using eqs. (30) and (31) of Brown, Koch and Rho [17|, one can convert 
our sigma-term attraction into an effective cr-exchange by 

V 4tt J eff 4tt 2pm K V ' 

where we have used Y*kn ~ 2m n and f = f n fa 93 MeV in the conversion. The agreement 
between the Miiller-Groeling et al value and ours is quite satisfactory. We thus see that 



4 T. Barnes, private communication. 
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there is evidence for attractive scalar interaction in the KN scattering consistent with 
Y, KN w 2m n . Q 

To see the role of the constants D and D 1 , let us first set them equal to zero. Then 
from (|17D and (|18|), we would predict (for Tikn ~ 2m u ) 



a|L^ ~ -0.30 fm, (24) 

KN ^KN — 3C ^nAce / c\ 

a I=0 = - -«0.46fm. (25) 

4vr/ 2 (l +m K /m B ) 

We have made clear in ( |25| ) that in the absence of the D and D' corrections, the af=Q 

is given mainly by the explicit chiral symmetry breaking. For this amplitude, the vector 

meson amplitude is zero. Equations (|24| ) and (|25| ) lead to ~ —0.11 fm which is too small 

compared with the empirical value (|2l| ) . Including the baryon-antibaryon pair contribution 

( |i~6| ) brings afj^ down to 0.40 fm, increasing ath to a~th ~ — 0.15fm. Now if in addition to the 

pair contributions, one decreases the value of T,kn by 30%, then one gets a ~ — 0.255 fm. 

Therefore it is not difficult to fit the amplitude a within the uncertainty of the parameters 

of the chiral Lagrangian. However it will then be difficult to understand in this way the 

scattering lengths aj=o anc ^ a f=i separately, even within the wide range of allowed values 

in n KN 
m a i=o- 

Suppose we ignore the D and D' terms but account for higher-order effects in the 
constants of the chiral Lagrangian. The most prominent quantity is /. Now at one loop 
(i.e., at v = 3), the kaon field will scale with the kaon decay constant fx, so / should 
be replaced by fx in the sigma term ~ Tjkn/ f 2 - On the other hand, the v = 1 term 
( |To|) involves vector- meson exchange and requires that / be identified with f w , the pion 
decay constant. Calculation of other v = 3 contributions has not yet been carried out, 
so it is difficult to quantify our argument but the use of fx instead of in the explicit 
symmetry breaking - which is equivalent to using (Exiv) e // ~ 0.69S^tv _ generally moves 
our theoretical values towards the central values of the empirical ones. In fact a somewhat 
lower value 

(^KN)eff ~ 0.6T, K N (26) 



would lead to an agreement between a t h and a, 



emp- 



5 The hypothetical higher-mass scalar repulsion introduced by Biittgen et al |8j in the K + N channel is 
probably simulating the so-called "counter-term" contribution to D and D' in the chiral Lagrangian discussed 
below and may be understood in the chiral bag language as a "van der Waals' repulsion" discussed by Vento 
et al Jl|]. 
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Let us now see what values of D and D are required to reproduce the "empirical" 



results (19). From (|7h and ([lj), we find 



D w 0.33/W - E KN /m 2 K , (27) 
£>' w 0.16/mjf - C/mlc. (28) 

For Y>kn ~ 2m7r, we have D ~ — 0.47/m_e with ms ~ 2rriK- The constant Z?' is ~ 0.46/m_e 
independently of the value of the sigma term. (If one takes afj* = 0, one gets D ss 
— 0.41/ms and D' 0.38/m# while for = —0.2 fm, one finds D ~ — 0.57/ms and 
D' ~ 0.54/ms.) Comparing with (|16|), one sees that the pair term accounts for only a small 
part of the constants D and -D', the latter disagreeing even in sign. The conclusion then is 
that there must be additional corrections to the scattering lengths at order v > 2. These 
corrections can be of two different classes: One class would be of the v = 2 terms coming 
from higher energy sector that is integrated out and is not saturated by the exchange of 
single resonances and the other would be loop corrections. As for the first class, we have 
no clue as to its mechanism. All we can do is to extract it from experiments. As for the 
second, since, as mentioned, there are no loop corrections to the D and D' terms, they must 
enter with v > 3, contributing effectively to the scattering lengths as do D and D' terms. 
Such calculations are in progress and will be reported elsewhere. 

That loops could be important either in -kN scattering [] or in KN scattering or in 
both is pretty much obvious. In irN scattering, our chiral Lagrangian (|5|) gives, at the v = 2 
order, the isoscalar ttN scattering length 

4n = \ ( a7T+P + aW+n ) = i^f\ l + m^/niB)}- 1 (2D^ml + S^) 

where S^^v is the ttN sigma term ~ 45 MeV and 

1 1 

D nN = -{dx+da) + ~(d 5 +d 6 ). 

Now if one takes the empirical value of at N [19|, al N = — 0.01m" 1 , then one must have 

Ar N » -0.27 fm. 

It should be noted, however, that the empirical value of a^ N can be reproduced by just taking 
the pair contribution alone (which is just a part of the D term in the Lagrangian), or, stated 
differently, by the pseudovector pole contribution fijl, or equivalently by the contribution 



6 This possibility was discussed in a somewhat different language by Delorme et al |J_9[ and also by K. 
Kubodera and H. Yabu (private communication). 
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of nucleon recoil pO. This means in our scheme that loop contributions must be present 



such as to substantially cancel the contribution of E^at to the scattering length IS]. This 
could be explained by a mechanism that enhances the pair contribution to ttN scattering 
by order ~ 2m n /m^ []. 

In what follows, we will simply use the empirical values (|27]) and ( pH ) effectively 
parametrizing v > 2 effects quadratic in kaon frequency. Later we will find that K~ con- 
densation is little affected by these v > 2 effects. The reason can be seen from (11). Since 
the D and D' multiply dtKdtK, these corrections scale with density as u^/rn^, where u> K 
is the kaon frequency. The loop corrections will also scale in the same way or what is 
more likely, even faster. The oj k decreases with increasing density, both because the kaon 
experiences the attractive scalar field of (|23|) and because of the attractive vector field from 
cj-meson exchange with the nucleon described by the first term of (|Toj) . Thus the terms 
prefixed by D and D' on the right-hand side of eq.(|ll|) (and loop corrections absorbed in 
them) decrease compared with the first term with the coefficient T,kn/ f 2 - Consequently 
although the influence of the Y*kn/ f 2 term on the scattering lengths may be cut by ~ 40 
% by the corrections, once the densities necessary for kaon condensation are reached, the 
influence of the correction terms will be greatly diminished. 

In relativistic heavy-ion collisions, in which there is not enough time for strangeness 
violation, KK must be created |22|. Interactions of kaons and nucleons via u;-meson ex- 
change have opposite signs for K and K, so that they average out. Interactions via p-meson 
exchange go out for isospin-symmetric matter. The baryon-antibaryon pair term and loop 
corrections do not decrease with increasing density compared with T,kn/ f 2 term, so that 
the corrections may become important. This issue would require a more careful treatment 
of higher chiral corrections. 



In-Medium Scaling 

So far we have been considering S-wave kaon interactions in free space. We would 

now like to take the chiral Lagrangian so defined and apply it to kaon interactions in nuclear 

matter. For this we have to account for the effect of the medium on chiral symmetry and 

other symmetries of QCD incorporated into the effective chiral Lagrangian. As argued by 

7 We conjecture that this enhancement may be understood by a mechanism analogous to the Nambu- 
Jona-Lasinio picture of the pion in terms of NN bubbles. In effective field theories such as ours, such a 
phenomenon could occur - as in BCS theory - through a collective mechanism which turns "irrelevant" 
terms (of higher chiral order) into "marginal" terms. See Polchinski ]2l]] on a discussion of this matter in 
the framework of effective field theories. 
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Brown and Rho H, the most economical way to implement chiral symmetry and trace 
anomaly of QCD - which are the most important properties of QCD at low energy - in 
nuclear medium is to endow to the in-medium effective chiral Lagrangian with the basic 
parameters of the theory that scale as a function of density. Arguments based on symmetries 
of QCD predict a universal scaling (valid at the mean field with the Lagrangian) of the 
quasiparticles relevant for the process in question S: 



m B rJ m *M rj f* 
m B ~ rn M ~ / 



(29) 



where the subscripts B and M stand, respectively, for baryons and mesons in the SU(2) 

sector (except for Goldstone bosons) and the star denotes in-medium quantities at a density 

p ^ 0. The unstarred quantities denote free-space quantities. At the mean-field level, the 

constant qa remains unsealed @. There is by now rather strong evidence that the scaling 

( p9| ) is valid in nuclei [24]. The consequence of this in- medium scaling on the effective 

Lagrangian (|l0|) and (|TT|) is simply that we replace / by /*. In medium, one is effectively 

including one-loop terms as mentioned above. 

As alluded above, to one-loop order in free space, the kaon decay constant fx receives 

0(Q 2 ) corrections relative to the pion decay constant f n . So working to that level in 

medium, fx should be distinguished from f n in the effective Lagrangian, not only in its 

magnitude but also in its behavior. This could be implemented with a Lagrangian consisting 

of Goldstone bosons and baryons only. One possible way is to write, following Gasser and 

Leutwyler |25|| , a symmetry breaking term that involves derivatives of 0(Md 2 ), where M 

is the quark mass matrix, and then perform a one-loop ChPT with this Lagrangian. This 

is very similar to what one does in skyrmion physics except that here only tree terms 

are considered. In this paper we shall not pursue this procedure. We find it far more 

transparent when the light vector mesons p, u, etc are explicitly present. Now in medium, 

we have f n — > /* « $(p)f n whereas model calculations (such as NJL) indicate that fx 

scales very little up to p w p , so we could assume f K w fx up to nuclear matter.^] This 

may be somewhat too naive but it should be more reasonable than taking it to scale like 

/*. The consequence of this argument is then that in medium / in eq.(10) is replaced 

by /* and that in eq. (pi]) by fx- The replacement in (1C) can be understood best in the 

description with explicit vector mesons. As noted in ||, the gauge coupling of the vector 

8 While there are no reliable lattice calculations in dense matter, we know what happens to the strange 
quark condensate (ss)o in high temperature. The calculation by Kogut et al j23| shows that the strange- 
quark condensate changes little as a function of temperature while the light-quark condensate does. This is 
consistent with our assumption. 
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meson g remains unsealed at the mean-field level, so the scaling of /* means the scaling 
of the vector meson mass m v by the KSRF relation m v = 2f* 2 g 2 . Thus the factor 1/f* 2 
just corresponds to 2g 2 /m v 2 for the vector propagator in the medium at zero momentum 
transfer. We will simply assume that neither T,kn nor D nor D' scales with density. In 
reality, there may be some scaling in both: the former because it is related to the condensate 
difference between the vacuum and the hadronic "bag" ; the latter because at least part of 
it may come from the pair term which depends inversely on the centroid baryon mass mg 
which presumably scales. We expect however that this is a fine-tuning that cannot be done 
accurately, so we will not pursue it any further. 

The remaining procedure for describing kaon condensation in neutron star matter 
is identical to what was done in [||, ||, so we will be very brief. Details are found in 
refs.[^, [|, ^6|. We simply mention that in medium, the Lagrangians ( JlOj ) and ( ]TT| ) generate 
one-loop terms and that kaon condensation will be triggered by the presence of electrons. 

Equation of State with Kaon Condensation 

The equation of state (EOS) describing the state of matter containing kaon conden- 
sates can be evaluated at the mean field level using the Lagrangian (||) with (|6|) and (j^), 
retaining all nonlinear meson interactions. For this, loop corrections may be quantitatively 
important. In contrast, evaluating the critical density is relatively simple as it is insensi- 
tive to nonlinearities and loop corrections. Here we give the results for critical density and 
composition of the condensed matter computed with a chiral Lagrangian consistent with 
the on-shell constraints to the chiral order v = 2. 

Before implementing the D and D' corrections, let us recall the previous results of 
refs.||, ^(J] obtained in the leading chiral order, v = 1 supplemented by a KN sigma term. 
The critical density is given as the nucleon density at which the pole of the kaon propagator 
D is equal to to the electron chemical potential in the absence of the condensate. The 
electron chemical potential is determined by the nuclear matter equation of state and the 
conditions of beta equilibrium and local charge neutrality. Equivalentlythe critical density 
is the nucleon density at which the energy density is lowered by the introduction of kaon 
condensate. In Table 1 [], we quote the results for the critical density without BR scaling 

9 With aim s — —67 MeV and a,2m s = 125 MeV determined by Gell-Mann-Okubo mass formulas, the 
parameter a?,rn s = —134, —222, —310 MeV corresponds respectively to Ekjv = Om^, 2.5m^, 3.8m^-. The 
case a^rris = —134 MeV corresponds to (ss)jv/(mm)jv ~ 0.07. This represents the least favorable condition 
for kaon condensation and will be used for illustrating the equation of state we predict. Note that the 
difference from eq.(^) is related to the well-known E^at problem in the tree order of chiral Lagrangians, 
which is resolved by going to higher orders in chiral perturbation theory. In the present problem, it is best 
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as given in refs.Q, |2(|, which are essentially those of Q but with the inclusion of muons in 
beta equilibrium with the neutron star matter. Here as well as in what follows, we use the 
potential contribution to the symmetry energy F(u) = u ||, 26, [27]]. 



TABLE 1 

The critical densities without D and D' terms and without BR scaling. The 
parameters held fixed are: a\m s = —67 MeV, a2m s = 125 MeV, fx w « 93 
MeV. 



F(u) a 3 m s (MeV) u, 





- 134 


4.11 


u 


- 222 


3.04 




- 310 


2.39 



Next we incorporate D and D' terms so as to be consistent (at the order v = 2) with 
the K + N scattering lengths. In what follows, we will simply take (|27|) and (28) valid to 
order v = 2. A more sophisticated analysis would require more accurate K + N data and 
loop corrections. The inverse propagator is 



D- l {uo) = uj 2 - m\ - n(w) 

(1 + [j 2 + (2* - l)^]up )u; 2 -ml + upj^-u + ^(£™ + (2x - l)C). 

(30) 



TABLE 2 

The critical densities in ChPT with v = 2 terms without (ti" s ) and with (uj! 1 ' 52 ) 
BR scaling. For parameter values, see Table 1. 



F(u) a 3 m s (MeV) 


u n c s 


< 


uf 


- 134 


4.20 


2.84 


3.15 


u - 222 


3.27 


2.49 


2.69 


- 310 


2.60 


2.16 


2.39 



to consider the numerical value of "E_r-jv" as a parameter. 
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The critical density computed with eq. ([?(]) without BR scaling is given in Table 1, denoted 
as u™ s . To implement the BR scaling in eq.(30), we replace f w — > /* and fx — > ~ /n- 
in the manner prescribed above. (We are disregarding the 0(Q 2 ) correction to fx in our 
numerical estimates.) As argued above, this corresponds to scaling the vector interaction 
only. Now if we take the scaling to be /* = /^/(l + 0.25u), we get the critical density 
labeled u^ 1 in Table 2. The results are only slightly changed if one uses a slower scaling, 
say, f* = f n /(l + 0.16n). They are given in Table 2 as u s2 . 

The composition of the condensed matter is illustrated in Tables 3 and 4. In Table 3 
are given the results obtained without BR scaling and in Table 4 those with BR scaling with 
/*//„-« (1 + 0.25u) _1 . We have shown only the results with the small T<kn ~ l.3m w (or 
o>z m s = —134 MeV) to illustrate the robustness of the prediction. The quantities reported 
in the Tables are: The chiral angle 9, the energy density gain Ae (in MeV/fm 3 ) , the 
chemical potential [i (in MeV), the proton fraction x, the kaon fraction xk = Pk/p, the 
electron fraction x e = p e /p and the muon fraction x^ = p^/p as a function of matter density 
u = p/po- 



TABLE 3 

"Nuclear star" composition without BR scaling for T,kn ~ Om^ {a^m s ~ —134 
MeV). 





u 




e 


Ae 


P 


X 


x K 


X e 




4 


20 





4 


0.0 


256.7 


0.195 


0.000 


0.111 


0.085 


4 


70 


31 


6 


-3.0 


223.9 


0.291 


0.180 


0.066 


0.046 


5 


20 


43 


6 


-11.2 


191.3 


0.360 


0.301 


0.037 


0.022 


5 


70 


51 


7 


-23.7 


159.9 


0.409 


0.380 


0.020 


0.009 


6 


20 


57 


4 


-39.3 


130.8 


0.444 


0.431 


0.010 


0.002 


6 


70 


61 


5 


-57.4 


104.5 


0.468 


0.464 


0.005 


0.000 


7 


20 


64 


5 


-77.1 


81.1 


0.486 


0.484 


0.002 


0.000 


7 


70 


66 


8 


-98.3 


60.4 


0.499 


0.499 


0.001 


0.000 


8 


20 


68 


6 


-120.5 


42.0 


0.509 


0.509 


0.000 


0.000 


8 


70 


70 





-143.5 


25.8 


0.517 


0.517 


0.000 


0.000 


9 


20 


71 


1 


-167.3 


11.3 


0.522 


0.522 


0.000 


0.000 


9 


70 


72 





-191.6 


-1.7 


0.526 


0.526 


0.000 


0.000 


10 


20 


72 


7 


-216.5 


-13.4 


0.530 


0.530 


0.000 


0.000 
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Some remarkable features in the results are as follows. 

The critical densities are robust with respect to the parameters of the effective La- 
grangian. The on-shell constraints bring only small modification. The in-medium 
scaling reduces the critical density from u c ~ 4 to u c ~ 3. 

With the in-medium scaling, the critical density is fairly insensitive to the numerical 
value of T,kn- 

Once kaon condensate sets in, the initially dense neutron matter turns quickly to 
nuclear matter, with the proton fraction x being already quite substantial just above 
the critical density This shows that with kaon condensates, a compact star is more 
likely a "nuclear star" rather than neutron star of the standard scenario. 

The kaon fraction becomes equal to the proton fraction slightly above the critical 
density; the proton charge is almost entirely balanced by the K~ charge. 

TABLE 4 

"Nuclear star" composition with BR scaling with f*/f w = (1 + 0.25u) _1 and 
^kn ~ l.3m w (a^m s « —134 MeV.) 



u 


e 


Ae 


A* 


X 


XK 


x e 




2.84 


0.0 


0.0 


209.6 


0.147 


0.000 


0.089 


0.058 


3.34 


28.0 


-7.0 


134.5 


0.356 


0.331 


0.020 


0.005 


3.84 


30.2 


-19.5 


92.8 


0.424 


0.419 


0.006 


0.000 


4.34 


29.9 


-33.3 


66.0 


0.457 


0.455 


0.002 


0.000 


4.84 


29.0 


-47.5 


47.6 


0.475 


0.474 


0.001 


0.000 


5.34 


27.8 


-61.7 


34.6 


0.485 


0.485 


0.000 


0.000 


5.84 


26.6 


-76.1 


25.3 


0.491 


0.491 


0.000 


0.000 


6.34 


25.5 


-90.5 


18.4 


0.495 


0.495 


0.000 


0.000 


6.84 


24.4 


-105.0 


13.2 


0.498 


0.498 


0.000 


0.000 


7.34 


23.3 


-119.6 


9.3 


0.499 


0.499 


0.000 


0.000 


7.84 


22.4 


-134.4 


6.3 


0.501 


0.501 


0.000 


0.000 


8.34 


21.5 


-149.3 


4.0 


0.501 


0.501 


0.000 


0.000 


8.84 


20.6 


-164.5 


2.2 


0.502 


0.502 


0.000 


0.000 


9.34 


19.9 


-179.9 


0.8 


0.502 


0.502 


0.000 


0.000 


9.84 


19.2 


-195.4 


-0.2 


0.502 


0.502 


0.000 


0.000 
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Conclusion: Astrophysical Consequences 

Our results show that the critical density is remarkably robust. It comes out to 
be u c ~ (3 — 4) quite independently of the constraints from kaon- nuclear interactions. 
What is also noteworthy is that once the BR scaling is implemented, the dependence on 
the strangeness content of the proton is considerably weaker. Even when one assumes 
that (JV]ss|iV) = which is certainly unreasonably conservative, the critical density is 
u c < 3. This robust nature of the condensation leads to some striking consequences on 
stellar collapse. 

Kaon condensation plays an important role in the collapse of large stars, once central 
densities exceed the critical density u c . Brown et al Q showed that electrons changed into 
K~ -mesons and neutrinos in a fraction of a second, the neutrinos being trapped for a longer 
time (of ~10 seconds). Since a large fraction of protons are present in kaon condensation, 
nuclear matter rather than neutron matter - formed in the conventional scenario - is reached 
and this has the effect of substantially softening the dense matter equation of state (EOS) 
because of lower symmetry energy of nuclear matter. Because of the greater binding energy, 
more energy will be emitted in neutrinos [26]. 



Effects of kaon condensation on the supernova explosion mechanism remain to be 
explored. However effects of kaon condensation on the structure of neutron stars and on 
some aspects of the explosion have been worked out by Thorsson et al pq] . 

Because kaon condensation softens the nuclear EOS at high densities, it substantially 



diminishes the maximum mass M max for neutron stars. Brown [28] found, with reasonable 



assumption about the compression modulus of nuclear matter that 

M max « 1.5M . (31) 

This has major implications for the formation of black holes in stellar collapse ||. Stars in 
the range of ~ 18 to 30 Mq can first explode, returning matter to the galaxy and then go 
into low-mass black holes of mass M BH ~ \.5M®. Upon collapse, stars heavier than ~ 30Mq 
drop directly into black holes without nucleosynthesis. These black holes are heavier, of 
mass M BH > WMq. Brown and Bethe || estimate that ~ 10 9 black holes have been formed 
in the above way in the galaxy. 
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Appendix 



Here we sketch briefly the derivation of the 1/uib corrections. Consider the lowest order 
Lagrangian with baryon fields, 



C = TiB (ip- m(l - i)B + F 7 ^7 5 K, B] + D 7 ^ 75 {yl M , B}) . 



(A.l) 



Rewriting the fields as B = -j^B a X a , V M = V£\ a and = A a ^\ a , where A a , a = 1,8 are 
the Gell-Mann matrices, the Lagrangian can be written 



C = B a ip ac B c - mB a 5 ac (l - i)B c + B a ^ 5 A a ^B c 



(A.2) 



Here 



ip 
A 



iS ac $ + i(2if abc )f 



ac 



F(2if abc )A b + D(2d abc )A b 



(A.3) 



and / , d are SU(3) structure constants. The equation of motion for the heavy baryon 
field is 

(g ac - m{\ - i)5 ac ) B c = (AA) 
with g ac = ip ac + 7^7 5 ^ c . Now, decompose B a using P + = ^(1 + i/) and P_ = ^(1 — if), 

B a = P + B a + P_B a = B% + Bt . (A.5) 



Applying P_ to eq.( A.4 ) from the left, we have 



Now applying P + to eq.(A/4) and using eq.( A.6 ), we get, modulo 0(1/ 



P+ ^ + 2^ 9aCP - gC ) B + = ° ' 
Thus the Lagrangian containing l/ms correction term is 

C = B%(g ae + ^-g ac P^)B e + . 

To write this explicitly, we note that the 0(A 2 ) terms are given by 

g acp_ g ce = ^^J^P-Yl^Af = -JFP+J? = ~(V ■ A aC )(v ■ A" 



(A.6) 



(A.7) 



(A. 



(A. 
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which leads to (m = tub) 



Prom the identities 

{2if abc ){2if cde ) 

{2if abc )(2d cde ) 
{2d abc )(2if cde ) 
{2d abc )(2d cde ) 

we have 

B a (v A ac )(v ■ A ce )B e : 



For S-wave scattering, 



Tr BA%B 



Tr(BA )Ti(A B) 



TrBBA 2 



TtBA BA 



2m b 



B a (v ■ A ac )(v ■ A ce )B e . 



(A.10) 



±Tr(\ a [\ b ,[\ d ,\ e ]]) , 
±Tr(\ a [\ b ,{\ d ,\ e }]) , 
±Tr(\ a {\ b ,[\ d ,\ e ]}) , 

^Tr (A a {A 6 , {X d , A 6 }}) - ^Tr(A a A 6 )Tr(A d A e ) , 



(A.ll) 



TrB (B 2 {v -A,{v A, B}} + F>F{v ■ A, [v ■ A, B}} 
+F 2 [v -A,[v A, B]} + FD[v ■ A, {v ■ A, 5}]) 

--B 2 Ti(Bv ■ A)Ti(v ■ AB) 
3 

(D + F) 2 Tr5(t> • A) 2 B + 2(D 2 - Y 2 )TiB{v ■ A)B(v ■ A) 



+(D - F) 2 TiBB(v ■ A) 2 - -D 2 Tr(Bv ■ A)Tr(v ■ AB) . 



(A.12) 



1 

4/ 
1 

4? 
1 

2/ 




^ {(NN){d t Kd t K) + (NfN) ■ {d t Kfd t K)) 



2 ({NN){d t Kd t K) + (NfN) ■ {d t Kfd t K)) 



— (NN)(d t Kd t K) 



The I/me corrections are therefore 



, 1 



£i = (NNdtKOtK) + D'--(NtN) • {d t Kfd t K) 



with 



/ 



Di = — 



1 



A8m B 
1 

48tob 



P 



((D + 3F) 2 + 9(D-F) 2 ) 
((D + 3F) 2 -3(D-F) 2 ) 



(A.13) 

(A.14) 

(A.15) 
(A.16) 
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